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In this article we replace the semi-heuristic derivation of the Rosenfeld functional of hard convex 
particles with the systematic calculation of Mayer clusters. It is shown that each cluster integral 
further decomposes into diagrams of intersection patterns that we classify by their loop number. 
This extends the virial expansion of the free-energy by an expansion in the loop order, with the 
Rosenfeld functional as its leading contribution. Rosenfeld's weight functions then follow from the 
derivation of the intersection probability by generalizing the equation of Blaschke, Santalo, and 
Chern. It is found that the 0-loop order can be derived exactly and reproduces the Rosenfeld 
functional. We further discuss the influence of particle dimensions, topologies, and geometries on 
the mathematical structure of the calculation. 
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I. INTRODUCTION 

Hard particle systems serve as reference fluids for soft, 
granular, and cellular matter. They interpolate the phase 
diagrams of molecular particles in the limits of low and 
high particle densities, where the influence of the smooth 
and attractive interactions is of secondary order. Phase 
transitions, as nematic and smectic, can be understood as 
entropic effects of the excluded volume 1 . This distin- 
guishes hard particle systems as the ideal starting point 
for perturbation theory. However, it also requires an ana- 
lytic representation, or at least detailed knowledge, of the 
free-energy functional. This requirement limits the use- 
fulness of computer simulations, as the minimization pro- 
cedure needs the functional form of the free-energy and 
not its function. To obtain a theoretical understanding of 
the liquid, crystalline, amorphous, and glassy states [2] 
we therefore need better analytical tools than are cur- 
rently available. 

During the last decades, several interesting approxi- 
mations have been developed to derive analytical expres- 
sions of the free-energy or the pair-correlation function 
[3]. However, most of them are restricted to hard spheres, 
such as the well known solution of Thiele and Wertheim 
[3HS] of the Ornstein-Zernicke equation in the Percus- 
Yevick approximation. A different approach was sug- 
gested by Reiss, Frisch, and Lebowitz [7], who used the 
result of Isihara and Kihara [5HTU] for the second virial 
coefficient for convex particles. Their scaled particle the- 
ory motivated Rosenfeld [TTHIS] to develop the funda- 
mental measure theory, which is based on the local de- 
coupling of the second virial integral, on the invariance of 
the free-energy functional under coordinate rescaling and 
on its solution of the scaled particle differential equation. 
The Rosenfeld functional is therefore a semi-heuristic re- 
sult, valid under the same assumptions as the scaled par- 
ticle theory. Nevertheless, its advantage is the explicit 



* stephan. kordenOrwth-aachen. del 



dependence on the particle geometry through the weight 
functions and its local representation of the free-energy 
functional as the convolute of weight densities. 

Later it has been shown by Rosenfeld and Tarazona 
P^^EO] that the functional leads to an inconsistency when 
the volume, filled with spheres, is restricted to a sin- 
gle layer, a tube or a one-particle cavern. This led to a 
geometrically motivated correction of the original form 
and resulted in a highly accurate functional for the fluid 
phase of hard spheres up to the freezing point [T71 |5T] . A 
different strategy made use of simulation data to go be- 
yond the Percus-Yevick approximation [531 US] ■ The sim- 
ple structure of the functional led to further applications 
for cylinders, discs, needles, and their mixtures f ^H - fSU] 
and alternative representations of the weight functions 
[5T1 [32] ■ For a recent review see also [53] . However, de- 
spite its success, it is not clear, how to extend Rosenfeld's 
approach further and how to go beyond the semi-heuristic 
construction of the functional. 

In this article we will begin an investigation to clarify 
the underlying mathematical and physical assumptions 
of the fundamental measure theory. In a first step it will 
be shown that the Rosenfeld functional is only the leading 
order of an infinite expansion of the free-energy in inter- 
section diagrams, which will be classified by their number 
of loops and intersection centers. The 0-loop order cor- 
responds to sets of particles that intersect in at least one 
common point and can freely rotate around this center. 
This intersection pattern corresponds to an infinite sub- 
set of Mayer clusters and will be derived in this work 
by generalizing the equation of Blaschke, Santalo, and 
Chern [5m38j . It will be shown that the infinite number 
of terms of the 0-loop contribution requires only the cal- 
culation of three Euler forms. The relation between the 
intersection probabilities and their corresponding subsets 
of Mayer clusters then allows the calculation of the free- 
energy functional via the virial expansion. However, in- 
stead of the virial series in the single-particle density, 
we have to interpret the expansion in Rosenfeld's weight 
densities. This reformulation of the virial expansion not 
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only reproduces the Rosenfeld functional as the 0-loop 
order but applies also to all further loop orders. 

The article is divided into two sections. The scope of 
part [n] is more general. Here, we introduce the concept 
of the loop expansion, |II A[ give some background infor- 
mation on differential and integral geometry, II B and 



lastly derive the weight functions from integral geometry 
in section [II C| Part |III| co nsiders the 0-loop contribution 
of the free-energy. |III A| recapitulates Rosenfeld's ideas 
leading to the semi-heuristic formulation of the funda- 
mental measure theory and Tarazona's corrections. This 
approach is compared in |IIIB| to our new ansatz via the 
virial expansion, where we derive the 0-loop contribution 
of the free-energy and prove its equivalence to Rosenfeld's 
fimctional. We end the article in IIVI with a discussion of 
the convergence of the loop expansion. 



II. THE INTERSECTION PROBABILITY OF 
PARTICLE STACKS 

A. Intersecting Particle Stacks 

So far, little attention has been paid to the approxima- 
tion scheme leading to Rosenfeld's free-energy functional. 
Instead, fundamental measure theory is based on a tower 
of three postulates that fix the functional's overall struc- 
ture: 1. the free-energy functional density is assumed 
to be a polynomial in weight densities, 2. uniquely de- 
termined by its the homogeneous scaling dimension and 
3. its solution of the empirical scaled particle differential 
equation. Actually, there is no physical argument that 
justifies these assumptions from first principles, and the 
failure of only one postulate could cause the downfall of 
the remaining parts. A first step in generalizing Rosen- 
feld's approach is therefore to test these three postulates 
and, if necessary, to replace them. This will be done 
in the following by comparing the third virial order of 
Rosenfeld's functional to its exact integral. 

Rosenfeld's truly remarkable step in developing a 
weighted density functional for spheres is the local split- 
ting of Mayer's f-function into weight functions and to 
recognize its relation to the Gauss-Bonnet equation and 
thus to the Gaussian curvature K [TTl [39l . For two 
particles Z?j , Dj , intersecting in a domain A = DiD Dj of 
coordinate vector rvi € Di n Dj , the f-function decouples 
into the convolute 
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depending on the particle positions fi £ Di, fj € Dj and 
the distance vector fij = fj — fi. In ([T]) the integration 
over the surface Sa has been transformed into a volume 
integral at normal vector n by equation ( |B12 1 and finally 
arranged in the symmetric form ([2|, assuming that the 
embedding space of the particles is of infinite volume. 

The derivation of the local decoupling of the f-function 
Q and its relationship to the Gauss-Bonnet equation will 
be explained in the following sections and deduced from 
the Blaschke, Santalo, Chern equation of integral geom- 
etry [34 37 . It provides an exact identity for the in- 
tersection probability of convex particles and determines 
the pre-factor l/(47r) uniquely. Based on this result, any 
Mayer diagram can be transferred into weight functions. 

The third virial cluster in this representation is now an 
integral over three intersection centers A, B, C, particle 
positions fi, fj, fk, and their corresponding rotations f2. 
Let us introduce the notation 
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for the differential volume element. The exact third virial 
integral has thus the form: 
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restricted by the loop constraint 

-TAB + rue + rcA = (5) 

of their distance vectors tab ~ t'b — ^a- Collecting 
terms according to their particle number and introducing 
Wertheim's 2-point density [4lti44] 
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for rAjTB G Di, equation Q can be written in the more 
symmetric form: 
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FIG. 1. The third virial Mayer diagram in the particle (left) 
and intersection representation (right): a) pairwise intersect- 
ing particles corresponding to the exact cluster integral and 
b) its approximation as the stack of third order. 



This is to be compared to the corresponding third virial 
integral obtained from Roscnfcld's functional: 
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which has a much simpler form, integrated over only one 
intersection center A <Z Di Dj C] Dk and depends on 
three weight functions only, instead of six in the exact 
expression Q. 

The principal difference between these two integrals is 
illustrated in FIG.[l] FIG.[I^) displays the generic inter- 
section pattern of the third virial integral with pairwise 
overlapping domains Q, whereas the corresponding fig- 
ure FIG. ^) shows the case of Q with only one such 
center. Rosenfeld's diagram is a degenerate third virial 
coefficient, obtained in the limit va = = fc: where 
the triangle of FIG. [I^) shrinks to the tree diagram of 
FIG. ^p)- The difference between the exact and the ap- 
proximated third virial integral is therefore the way in 
which the particles intersect each other. 

Instead of the graphical representation of intersecting 
particle domains, it is sufficient to symbolize the inter- 
section patterns in "intersection diagrams", where the 
particles correspond to lines and intersection centers to 
the position where the lines join. The corresponding dia- 
grams of the third virial are shown on the right of FIG. [T] 

Rosenfeld's functional contains an infinite number of 
further virial contributions. These are obtained by Tay- 
lor expanding the singular parts in powers of the weight 
function ujy and have the generic form: 
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corresponding to Mayer diagrams, whose intersection do- 
mains have been contracted into one single domain, as 
shown in FIG. [2] However, only completely connected 
Mayer clusters interact in such a way that each parti- 
cle interacts with each other. Rosenfeld's functional is 
therefore the sum over an infinite number of completely 



FIG. 2. The Rosenfeld functional is the 0-loop approximation 
of the free-energy. Each intersection diagram corresponds to a 
completely connected Mayer cluster, contracted into a stack. 
The sum over all such diagrams is symbolized by a crossed 
circle. 



connected diagrams that are further contracted into one 
intersection point. 

The arguments, obtained so far, can be summarized 
in the following way: The exact free-energy functional is 
not representable by Rosenfeld's weight densities alone. 
Instead, the third virial integral ([t]) is a function of 
Wertheim's 2-point densities ([6]), and it is natural to as- 
sume that this result has to be generalized to arbitrary 
k-point densities. Next, as the Mayer function ^ is itself 
invariant under coordinate scaling, it is not possible to re- 
strict the functional form by its scaling dimension. From 
this follows that the three postulates of FMT, including 
the empirical scaled particle differential equation, have 
no deeper physical basis. On the other hand, we have 
also seen that Rosenfeld's functional approximates and 
re-summes a certain class of Mayer diagrams contracted 
to one intersection point, as shown in FIG. [l] This of- 
fers an alternative approach to derive the functional and, 
most importantly, it also opens a path to derive higher 
order corrections. 

The central object of FMT is the sum of contracted 
intersection diagrams, shown in FIG. [2] Because of its 
importance, let us introduce the name "stack" for indi- 
vidual parts and "universal stack" for its sum, defined 
by: 

Definition 1 A stack of order k — ord{Stk) is a set of 
i — 1, . . . , fc domains Di, intersecting in at least one com- 
mon point and free to translate and rotate around this 
center: 
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The universal stack is the formal sum over all stacks in- 
tersecting at the same point 
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In the following sections we will prove that the inter- 
section probability of the universal stack $0 reproduces 
Rosenfeld's functional 
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However, this is only the first hint to a more general 
structure: When completely contracted intersection dia- 
grams correspond to a free-energy functional at low pack- 
ing fraction, it is natural to assume that diagrams, not 
completely contracted, provide higher order corrections. 
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FIG. 3. Mayer clusters of the fourth virial order, translated 
into intersection diagrams and ordered by the tuple {g,n): a) 
(3,6), (2,4), (0,1), b) (2,5), (1,3), and c) (1,4). 
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FIG. 4. The displacement of one particle causes a shift in 
the position of all neighbors that are in direct and indirect 
contact, resulting in long range correlations between particles. 



Rosenfeld's functional is exact for the second virial or- 
der. The third virial integral, however, is only an ap- 
proximation, as shown in FIG. [T] Adding the exact third 
virial diagram will therefore result in an improved func- 
tional, corresponding to three additional intersection cen- 
ters and the loop constraint ([S]). In principle, it is possible 
to add arbitrary intersection diagrams to the functional, 
systematically derived from the Mayer clusters. As an 
example consider FIG. |3] where the 4-particle Mayer dia- 
grams are shown together with their corresponding set of 
intersection diagrams and contractions, ordered by their 
number of loops g and intersection centers n. This clas- 
sification by the tuple {g, n) comes natural as the calcu- 
lational complexity increases with both. However, they 
have also a direct physical interpretation. 

The loop order g counts the number of constraints, 
restricting the coordinates of intersection domains: 



^41^2+- 

rB,B2+- 



rBi^^Bi = 



g loop constrains (13) 



for g loops with k,l, . . . intersection centers. In this way, 
correlations are generated between particles that other- 
wise do not interact directly via a potential function. 
This distinguishes the zero loop order g — where no 
such constraints exists, providing a plausible argument 
why Rosenfeld's free-energy functional describes only the 
fluid regime below the first phase transition and predicts 
a maximum packing fraction at = 1, independent of 
the particle geometry. The solid phase region, on the 
other hand, requires long range correlations between par- 
ticles, such that shifting one particle leads to the displace- 
ment of others, as shown in FIG. |4] 

These considerations make the number of loops g and 
intersections n convenient indices to group the diagrams 
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FIG. 5. The re-summed and regularized third virial integral; 
Adding additional particles to an existing intersection point 
does not increase the calculational effort. 
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FIG. 6. Going beyond Rosenfeld's functional: A first approx- 
imation for the 1-loop order contains the re-summed second 
virial integral and the regularized third virial integral. 



and to define the "loop expansion" of the free-energy ex- 
cess functional: 
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where each element corresponds itself to an infinite 
number of intersection diagrams. Examples are shown for 
$0,1 = $0 in FIG. [2] and for $1,3 in FIG. [Sj It is worth 
pointing out that some of the contracted 4-particle dia- 
grams of FIG. [3] turn up as corrections of the second and 
third virial order. Actually, it will be shown in the follow- 
ing sections that the calculational effort does not increase 
when additional particles are added to an existing inter- 
section point. Any individual intersection diagram can 
therefore be replaced by a "resummed diagram" , with 
each intersection point replaced by the universal stack. 
Resummation is therefore a central aspect of FMT, as it 
generates the pole structure in the free volume 1 — n^, 
which is so characteristic for Rosenfeld's functional. 

A natural extension of the current functional is the 
combination = *&o.i + ^1.3, shown in FIG. [6] 

However, as parts of $0,1 are already included in $1.3, it 
is necessary to "regularize" the loop diagram by exclud- 
ing the case Itab] = l^scl — \^ca\ — 0, shown in FIG.jl] 
where the distances between the intersection coordinates 
vanish. All loop diagrams are understood in this way, 
excluding the case of collapsing loops and thus ensuring 
that regularized diagrams are uniquely defined. 

Apart from the resummation of intersection points, it 
is also possible to sum up diagrams of identical loop or- 
der. One example is displayed in FIG. [7] The analytical 
structure of the generating function $1 can be derived 
from the virial expansion, as the ring diagrams are for- 
mally identical to Mayer clusters. With the symmetry 
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FIG. 7. The unregularized free-energy at 1-loop order is the 
sum over the second virial and all further ring diagrams. 



factor {k — l)!/2 for a ring of k particles and the sim- 
plifying notation C^^uja(jJb for the f- function the 
1-loop free-energy yields the formal expression: 
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where the angular brackets indicate the integration over 
the coordinates. The 1-loop free-energy contribution is 
therefore of a completely different structure than Rosen- 
feld's functional, signaling a logarithmic divergence, de- 
pending on the particles' geometry. 

Having identified the approximation scheme behind 
Rosenfeld's functional, we will now begin with the de- 
velopment of the mathematical framework necessary to 
derive the intersection probability of the universal stack. 



In this way, the hypothesis ( 12 ) will be proven by direct 
calculation, which is the basis for the resummation of in- 
tersection points and all further constructions that will 
be considered in following papers. 



B. Some Relevant Information on Differential 
Geometry 

1. Intrinsic Geometry 

The derivation of the intersection probabilities requires 
the introduction of some mathematical conventions [451 - 
ITT] and the discussion of physical constraints. 

Let D denote an Euclidean, Riemannian manifold of 
3 dimensions, sufficiently differentiable to allow for the 
calculation of the Euler form. Manifolds of this type in- 
clude a variety of geometries as convex and concave par- 
ticles, Klein's bottle, tori, polyhedrons, cylinders, hollow 
spheres but also non-compact structures. The mathe- 
matical requirements are therefore not very restrictive. 
However, we also have to take into account the physical 
constraints. In the formulation of Mayer's f-functions, 
the cluster integrals determine the intersection probabil- 
ity between particles. However, the physical particle do- 
main is not only restricted by its surface. Instead one 
has to determine the region that is inaccessible for other 
particles. FIG. [8] shows two examples, where the cor- 
responding mathematical intersection probability is zero 
but not its physical one. FIG. [8^) shows two linked tori. 
In a fluid of single tori, such a configuration has to be 
excluded as the particles cannot penetrate each other. 
The same applies to the system of a particle inside a 
concave domain, whose opening is smaller than the par- 
ticle's smallest diameter, as seen in FIG.[8}d). When the 
geometry of the first does not allow to enter the inner 
region of the second, it has to be excluded, i.e. counted 
as part of the domain of the latter particle. 

Although physically related, the mathematical nature 
of these two examples is very different. The case of two 




FIG. 8. Physical constraints restrict the regions accessible for 
other particles. Forbidden configurations are: a) linked tori 
and b) particles inside a cavern (bottle) with pore opening 
smaller than the particles' diameter. 



tori is related to Euler's linking number 133] and be- 
longs to the topological class of honiotopically non-trivial 
intersections. Another example is the intersection of hol- 
low spheres as realized by fuUerenes. Both cases are 
related to topological classes that follow by successive 
variations of Euler forms [SD]. And although they are 
not required for the current article, they are interesting 
enough to give a short account further below. The sec- 
ond case FIG. [8}d) is more difficult to solve. Here, we 
have to introduce a fictitious membrane at the opening 
of the pore, whose surface vector is always antiparallel 
to the surface vector of the docking particle. Such con- 
figurations lead to the vanishing of certain contributions 
of the intersection probability between particles, as we 
will show in the next section, and might give new insight 
into the isotropic-nematic phase transition. Because of 
these additional complications, we will exclude homotopi- 
cally non-trivial particles as well as concave geometries. 
The discussion simplifies further, when boundaries are 
excluded, leaving us with 3-dimensional convex particles 
embedded into the flat Euclidean space M"^. 

The geometry of a physical particle depends on in- 
trinsic and extrinsic properties, i.e. the properties in- 
dependent and dependent on the embedding. It would 
be therefore sufficient to consider 2-dimensional surfaces 
and their embedding into IR'^. However, at this point it is 
worthwhile to discuss Cartan's formulation of differential 
geometry [45l [46l [51] for general dimension, as some of 
the results degenerate for low dimensional spaces. 

Let the particle S be a rt dimensional, orientable, dif- 
ferential Riemannian manifold without boundary. Sup- 
pose further that the manifold can be covered by a set of 
open coordinate patches S = UaUa, each one isomorphic 
to ]R" and labeled by a local, orthonormal coordinate 
frame (p, e^""*, . . . , el""*) at the point p G L/q. The lo- 
cal frames at overlapping regions Ua n Ujs are related to 
each other by differentiable coordinate transformations 
gapip) ■ Ua Ci Up ^ SO(n). The matrix valued tran- 
sition functions are invertible = gpa and fulfill 
the cyclic condition ga/^g/^-yg-ya = 1 at triple intersections 
Ua D Ui3 O U-f. These preliminaries define the tangential 
bundle TE with the local section {p,ei, . . . , e„) S F(TI]) 
and the cotangential bundle T*Y, as its dual space, re- 
lated to TE by the metric of M" 

ei ■ ej = i-jij (16) 

and its differential structure. The vielbein 0, and con- 
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nection forms are defined by 



dp 



dci = uJijCj for p G E 



(17) 



and transform under the coordinate change e[ — gijCj as 



= Ski^kigi] + 9kt^dgkj 



(18) 



where summation over paired indices is understood. The 
connection is therefore not a tensor and can be locally 
replaced by a trivial gauge. 



The vanishing of the second exterior derivative of (17) 
defines the torsion and the curvature form 

T, = dOi - A 6j , Qij = dwij - LOik A Wfej , (19) 

which transform as a first and second rank tensor 



(20) 



The constraint = of a Riemannian manifold is there- 
fore independent of the coordinate system and introduces 
a global relationship between the vielbein and the con- 
nection forms. 

Torsion and curvature carry local information about 
the geometry of a manifold, always restricted to single 
coordinate patches Ua and depending on the chosen co- 
ordinate system. Globally defined forms, on the other 
hand, are necessarily invariant under coordinate trans- 
formations. An important class of such functions was 
introduced by Chern [171 in extending the notion 
of class functions f{g~^xg) = f{x) from group theory. 
From ( 19 1 follows that the curvature form transforms un- 
der the adjoint representation of SO(n). Natural choices 
are therefore the determinant and the trace of a polyno- 
mial in f2, whose differential form is of the same order as 
the volume form of S or any submanifold thereof. Chern 
defined the Euler form or Euler class 
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for even dimensional manifolds n 
as the Euler characteristic 



2k and its integral 
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with the normalization chosen such that its result is 
whole- numbered for the sphere x{S'^) = 2 and the g- 
holed torus x{Tg) = 2 — 2g. The integral is a topological 
invariant and central for many areas of mathematics and 
physics |49j . It is therefore not surprising to discover that 
the Euler form also enters the discussion of hard particle 
physics as the intersection probability of particle stacks. 

The Euler class is the highest possible form for even di- 
mensional manifolds from which derives a series of invari- 
ant differential forms by successive variation 5ui = g~^dg. 
The resulting Chern-Simons classes |47l[5Q] determine the 
failure of the form to be invariant under the coordinate 



transformations (18). As an example consider the case 



of n = 2 dimensional manifolds with the transition func- 
tion g = exp (iA) G U(l). The curvature reduces to the 
exterior derivative ^ — duj and its variation 5u) — idX to 
a iM. valued function: 
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When the first integrant is rewritten by the Gaussian 
curvature if, the second by the geodesic curvature Kg and 
the last integrant by the interior angles, we obtain from 
( 22 1 the Gauss-Bonnet equation for the 2-dimensional 
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with non-contractible curves along 91] and additional 
vertices at the singular points. To get a better under- 
standing of the origin of these additional contributions, 
remember that the Euler form counts the angular change 
of the normal vector, while moving over the surface of the 
embedded manifold. For smooth, Riemannian surfaces 
this is always An, but boundaries and singular points 
contribute additional angular changes and generate the 
Chern-Simons terms. 

It can be shown [371[SD] that the two equations fi = dw, 
6u! = idX generalize for the Euler form Qjfc •= Pf(f^) for 
arbitrary even dimension to a sequence of characteristic 
classes 



rAr\n = (25) 



for m = 1, ... ,2k. Each variation now produces a new 
characteristic form of one order less than its predeces- 
sor. And in the same way as the geodesic curvature is 
an invariant form for the 1-dimensional curve 9E, it is 
natural to apply the odd differential forms of Q^fe-m-i 
to odd dimensional manifolds of non-trivial homotopy 
group. Euler's linking number and the intersection num- 
ber of hollow spheres are special cases of these forms. 
In the notation of 50|, they correspond to Q2 and Q3 
and derive from the Euler class of a 4 and 6 dimensional 
manifold. However, for convex geometries, which we will 
consider in the following, it is not necessary to take these 
classes into account. 

Apart from the geometric interpretation of a Rieman- 
nian manifold, there is also the relation to Lie groups, 
whose vielbein and connection forms constitute the basis 
of a Lie algebra [351 SSI 123] represented by the matrix 



aA = \ Q ] e iso(n) 



(26) 



whose elements satisfy the Maurer-Cartan equations 
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They are related by the inner derivation isCTA = 
aAiXs) = TjAB to the more commonly used commuta- 
tion relation [Xa,Xb] = Cab-^c and Jacobi identity. 
The corresponding Lie group is the Euclidean or isomet- 
ric group that locally splits into the semi-direct prod- 
uct ISO(n) = SO(n) k E„ of rotations and translations. 
Its Lie algebra elements ujij and 9i transform under the 
mapping ( |18[ ) and span a n{n + l)/2 dimensional space 
consisting of the n{n — l)/2 connection and n vielbein 
forms. 

The integral over all rotations and translations is there- 
fore related to Haar's measure of the isometric group 



= /\'i=i^in A-tTi a;j„_i A . . . A W12 A dvol(£'„) 
= dS""-^ A dS''-^ A...AdS^ A dvol{E^) 
^ dvol(SO(n)) A dyol{En) , 

where we made use of the coset representation: 

SO(fc)/SO(/c-l) = S''=-i . 



(28) 



(29) 



Evaluating the integral yields then a product of volumes 
of spheres, with values: 



Ok = voliS' 



27r^ 

rC^) 



(30) 



whose first elements are Oi — 27r, O2 = 47r, 



2. Extrinsic Geometry 

Up to now, we have only considered the intrinsic prop- 
erties of the particles' geometry. However, the move- 
ment in a background space requires the choice of a suit- 
able embedding. For physical reasons it is natural to 
consider the flat Euclidean space and to imbed the k 
dimensional particle, e.g., into the first coordinate di- 
rections of the local frame (ei, . . . , efe, e^+i, . . . e„) with 
the corresponding nontrivial coordinate transformations 
ISO(n)/(ISO(fc) X ISO(n - fc)). To avoid the additional 
problems that occur when discussing this complicated 
coset structure, we will restrict the dimension of the em- 
bedding to fc = n — 1. The group ISO(l) consists then 
only of the translation in one direction and can be ex- 
plicitly separated in the following equations. 

This choice is the simplest possible embedding and at 
the same time also the physically most relevant one. The 
manifold D ^ M" is now a n-dimensional domain in E" 
and bounded by its surface dD. Following the outline 
of [36], we choose the outward normal direction of the 
surface to point along e„, such that the tangential direc- 
tions of dD correspond to the first n — 1 elements of the 
local frame (ei, . . . , e„_i, e„) of M". The corresponding 
directions are differenced by the index convention: 



i,j = l,...n and a, /3 = 1, . . . , n — 1 . 



(31) 



The associated Pfaff system ^51 [5T] of the integrable 
submanifold is then defined by the constraint 



: on dD 



(32) 



Applied to the vanishing torsion of the Riemannian mani- 
fold 



= d0„ 



hafsOp A I 



KaVa A Va , (33) 



it allows an algebraical solution of the equation by the 
symmetric matrix /iq^ and to define the principal curva- 
tures Kq, and principal vectors as its eigenvalues and 
orthonormal eigenvectors. In the form, de„ = HaMae-a, it 
is also known as Rodrigues formula. 

Splitting the 71-dimensional curvature il^"^ into normal 
and tangential directions 



dl^r, 



■ LOaP A UJj3ri 



(34) 



(n-l) 
af} 



A UJ< 



yields the Gauss and Gauss-Codazzi equations '45] , fur- 
ther reducing to 



(n-l) 



and 







(35) 



in the case of fiat embedding. The first equation relates 
the intrinsic curvature of the particle to the normal con- 
nection forms of the embedding. Whereas the vanishing 
of ujna under the n — 1 dimensional covariant derivative 
ensures the decoupling of the normal coordinate trans- 
formations from the tangential ones; the forms 

are 

therefore horizontal |45| . without the need of introducing 
equivariant differential forms |54| . 

In the definition of the embedding we have assumed 
that the normal vector e„ points outward from the com- 
pact particle surface. This corresponds to a special 
gauge choice in the 0{n) coordinate transformations of 
M" and restricts the group to SO(n). But this local 
gauge does not extend globally, where both orientations 
Z2 = 0(n)/SO(n) have to be taken into account. The 
Euler characteristic, derived by the intrinsic curvature 
([19]) and by Gauss's equation (35 1, wiU correspondingly 
differ by a factor of two 



x{dD) = 2x{D 



(36) 



The kinematic measure ( 28 ) of an embedded particle 
of odd dimension n = 2k + 1 can now be calculated by 
combining (21 [22 35 36) and observing that the nor- 
malization of the Euler characteristic is p roportional to 
02k = 2(47r)'=A:!/(2fc)!, as follows from ^ 



. 2/c+l . 2fc+l 



= J Pf(rj) A (ivol(SO(2fc)) A dvoliE2k+i) (37) 
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= x{dD) vol(SO(2fc + 1)) vol(£;2fe+i) . 

For a 3-dimcnsional manifold in M"^, the corresponding 
integral reduces to 



U)l2U!i3U)23Ol0203 = / KiK2i^lJ^2- / 0112 (ivol(M^ ) 



J K{dD)dS A -271 V 



K{dD)d{nfA)(frA 



(38) 



with the volumes of vol{M.^) = V and vol(SO(3)) = 8tt^. 

Note that the kinematic measure of a Riemannian 
manifold would vanish for dimensional reasons, as the 
vielbein and connection forms are not independent. It is 
therefore necessary first to interpret the integrant as the 
Haar's measure and only afterwards to incorporate the 
geometric constraints. 

This equation is of course closely related to Chern's 
original derivation of the Euler class [55] ■ Here however, 
the difference lies in the relation between geometry and 
isometric group, which focuses on the alternative inter- 
pretation as the kinematic measure of a particle, moving 
in a flat background. For two intersecting particles it thus 
determines the intersection probability, averaged over all 
rotations and translations. It is therefore identical to the 
second virial integral and explains the appearance of the 
Gauss-Bonnet equation p4| ) in the calculations of Isihara 
and Kihara ^ , Rosenfeld |39j , and Wertheim [?T] . 



C. The One, Two, and Three Particle Intersections 



1. Comments on Integral Geometry 



The generalization of (37) to two and more intersect- 



ing particles leads us into the field of integral geome- 
try, whose differential geometric formulation goes back 
to Minkowski Weyl Blaschke |34|, Santalo [35], 
and Chern [55H55] , who observed that the invariant forms 
of integral geometry can be traced back to the Euler class 
( 37 1 . One intriguing result is the fundamental kinematic 
equation j35j 



Vn 

In 



j x{Di n gD2)dg = 



\ ^ / n 



ISO(n) 



^ \k 

fc=0 



Mfe(A)M„-fc(^2) 



1 



7„ = VOl(S0(7l)) , Vn = -On 



(39) 



and the observation that the coupled geometry of two 
intersecting manifolds reduces to a simple pairwise prod- 
uct of Minkowski measures or integrals of mean curvature 
Mk- For n = 3 it reproduces the equation of Isihara and 



Kihara of the second virial coefBcient. Actually, they 
used for their calculation an early result of Minkowski 
[55) . In fact, it was the starting point for our current 
investigation and offers a direct, albeit less general, ap- 
proach of deriving the intersecting probability, which is 
why we have added their calculation in a somewhat clar- 
ified form in appendix \K\ 

There are several ways to derive the fundamental kine- 
matic equation ([39]). Probably the simplest one uses 
the expansion of the Steiner polynomial |35| . another 
one Blaschke's cut and paste construction [31| of sub- 
spaces. The most fundamental, although more elaborate 
approach is Chern's explicit derivation |36| of the Eu- 
ler class from the kinematic measure (37). Its advan- 



tage is the explicit local formulation in connection forms 
that will be important for its decoupling into Rosenfcld's 
weight functions. This ansatz is therefore the natural 
starting point for relating Rosenfeld's approach to inte- 
gral geometry. 

The generalization of (37) to a particle stack Stfc+i is 



easily achieved but requires some normalization to get a 
well defined result. First, we have to fix the position and 
orientation of one particle in St^+i to remove the volume 
dependence on the embedding space V = vol(M'^), gen- 
erated by moving the stack in the background manifold. 
Furthermore, it is useful to define the kinematic measures 
of the particle domain D and its surface S 



dD = a;i2a;i3a;23 A 6'i6'203 
dE = wi2 A 6*1 6*2 



(40) 



analogously to ([3| . The kinematic measure of ( 37 ) or ( 38 ) 
generalizes then to the integral average of A: -|- 1 particles 



1 1 



in 



dDi A ... A dDk+i 



Stk+1 



K{dStk+i)5{nrA)d''rA 



(41) 



X ^72 A . . . A d7fc 



-1-1 : 



with the Gaussian curvature K integrated over the do- 
main A = Stfc+i at fixed kinematic measure 7, as defined 
in (|3|, and integrated over the center of gravity, repre- 
sented by 7i. 

The boundary of the stack 9Stfc+i can be determined 
by the algebraic relations of the homology operator [15] . 
As an example, consider two intersecting manifolds that 
itself have no boundary d^D — 0. The application of d 
to the second order stack 

d{Di nD2) = dDinD2 + Din dD2 + dDi n dD2 . (42) 

is thus a sum of intersections, wherein each successive 
application of d reduces the dimension by one. This re- 
stricts the possible number of boundary operations to 
the dimension of the embedding space by the constraint 
d'^X = for any 3-dimensional manifold X. The infinite 
number of virial contributions, shown in FIG. [2] reduces 
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therefore to the derivation of three Euler forms, corre- 
sponding to one, two, and three particles. 

The calculation of (41 1 can be further simplified by in- 



cluding the physical constraint of indistinguishable par- 
ticles. To obtain the correct combinatorial pre-factors, 
let us define the formal sum 



M 



(43) 



of 1-particle densities and domains. It is the homologous 
operator of Rosenfeld's weight densities and parallels the 
notion of a divisor in algebraic geometry. The repre- 
sentation of the free-energy functional in n reduces the 
problem of determining the boundary of the stack of /c-f 1 
different particles 9Stfc+i to the corresponding analysis 
of a stack of fc + 1 identical manifolds, whose boundary 
reduces to a sum of three terms 

astfc+i = (fc + 1) E n stfe + (/c + i)fc E n s n stfc_i 
+ (fc + i)fc(fc -i)EnEnEn Stfe_2 

(44) 

in the shorthand notation S = dD. Using the linearity 
of the Euler form and its vanishing for odd dimensional 
manifolds, it translates to the corresponding Gaussian 
curvature 



K{d'&tk+i) = {k + i) k{j:) + (fc + i)fc k{j: n e) 
+ (fc + i)fc(fc - 1) K{T, n E n E) 



(45) 



that will be derived in the following. The first two terms 
are known from Chern 36 , who obtained the result for 
two intersecting manifolds of arbitrary dimension. An 
independent approach was used by Wertheim [41] . How- 
ever, the three-particle intersection is new and will be 
presented parallel to the summary of the previous two 
cases. The corresponding generalization of Chern's ap- 
proach to an arbitrary number of particles and dimen- 
sions has been developed in [S^ and will now be applied 
to three dimensions. 



2. The One Particle Euler Form 

Let us begin with the simplest case KCS) of one parti- 
cle, moving in a background of fc domains. Following the 



derivation of ( 38 ) the product of the connection forms can 
;he principal basis a;i3AcLi23 



be rewritten in t 
reducing the kinematic measure of E 



KlK2i^lAl/2, 



1 I f , 11 

— - — / dDi = — ^ — 



'^13^^23'^12P1P2C'3 



r(_Di 



47r 



A=Dx 



(46) 



with the Gaussian curvature kg and a factor of 27r from 



the integral over a;i2. The first part of the integral (44) 



for a stack can now be written as 

^KGidDuTA) 6{nrA) e(i?2, r4) • • ■ 6(^^+1, r^) 



X d ^72 . . . d'jk+i 
oj'^^^oji^K..JJ'+^^dhAdj2.-.djk+i 



(47) 



in the weight functions 



OJ. 



1 

in 



KGSidD,), 4) = e(A), (48) 



where the integration domain has been formally extended 
to the complete embedding space V by the Dirac- and 
Heaviside- function S and 8, with 5{dDi) understood as 
restricting the volume integration to the surface 5{nrA) 
at the intersection point va G St^:. 



3. The Two Particle Euler Form 

Some more efforts requires the derivation of the sec- 
ond Euler form if (Ei n E2) that determines the angu- 
lar change between the two normal vectors at the 1- 
dimensional intersection submanifold. It parallels the 



geodesic curvature Kg of the Gauss-Bonnet formula ( 24 ) 
and can be seen as the real space generalization of the 
Chern-Simons class. Its derivation begins with the con- 
struction of a proper coordinate system at the intersec- 
tion space. Let us introduce the bases Ei : (61,62^^63^'') 
and E2 : (61,62 ,63 ) with the common direction Ei n 
E2 : 61 along the 1-dimensional submanifold and the in- 
tersection angle 

cos (012) = (e^^^e^^)) for < 0i2 < 27r . (49) 
Following |57j . we define the intersection determinant 

M,=det(6«6i^))|,t,.^i (50) 
for fc intersecting surfaces. The first two cases are: 



M2 = 1 - 612 = §12 

Af3 = 1 - C12 - Ci3 - C23 -I- 2C12C13C23 , 

where we used the shorthand notation: 



sin (0ij) , dj := cos {(t>ij) 
5(7) := sin (7) , 0(7) := cos (7) . 



(51) 



(52) 



The local frame of the intersection manifold in W is 
spanned by the vector field 61, 63^'*, 63'^'' for 0i2 7^ 0, from 
which one obtains an orthonormal basis by the Gram- 
Schmidt process 

(1) 
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«2 
Vi 



1 

Sl2 

ei . 



(53) 



As explained before, the Euler characteristic counts the 
angular change of the normal vector, while moving from 



Cg^-* to e'^' . To interpolate between those two vectors, we 
introduce a SO (2) rotation in the range < 7 < (/)i2 



= (2) 



% = c(7)w3 + 5(7)^2 

m = ~s{j)v3 + C(7)W2 

771 = ei . 



(54) 



One of the two equivalent vectors, 773 or 772, is now the 
new outward pointing normal direction. Let us chose 
773 and derive the corresponding Euler density for the 
intersection Si n S2: 



mdva A ri2dri3 
_ 1 
S12 



S(012 - IM3 + 5(7)0; 



^13 



A 



(55) 



with the definition = eide^^ of the new connection 
forms for the particles i = 1,2. Integrating over 7 



771^773 A 772^7^3 = 



1 - C12 



,(1) 



■^13 \ 



S12 - " (56) 

= if (El n i]2)rf(Si n E2) 

yields the differential Euler form. Observe, that the an- 
gular dependent factor (1 — Ci2)/si2 remains finite even 
in the limit of anti-parallel vectors 4>i2 — >■ tt when the 
remaining kinematic measure is included, which will be 
derived below. 

At the intersection Ei n E2, a S0(2) transformation 



(18 1 relates the vector frames of the two particles 



= (2) 
'3 

(2) 
--2 

.(2) 



21263^^ 4 

-^1263^' 



(2) 
51262 

+ ci2e^^^ 



and the boundary condition f?3*''|aEi — their corre- 
sponding Pfafhan systems ( 32 ) . The transformed differ- 
ential forms 



(57) 



UJ 



g(2) 

(2) 
13 

,(2) 



Cl2^^U 



'"23 



512^'^ 
- Si2^'l^2 



)12 



(58) 



are therefore understood modulo 9'^\uj^^l. With these 
relations, the reduced kinematic measure of Ei n E2 can 
be derived, with the first particle fixed in the embedding 
space and the second one free to move: 

d(Ei n E2) A dD2 

- A fl(2)fl(2)/l(2), ,(2), ,(2) (2) 

— e^i A fi ^2 L0^2 LO^^ LO23 



2 /\ p(2)fl(2),.,(2) 



^12 



(512) 

(si2)^ dEi A dE2 A I 



■^12 



A d(t>i 



(59) 



^12 



with the kinematic measure of the surface defined in ( 40 1 



The decoupling of the Euler form ( 56 ) and the kinematic 



measure ( 59 1 for two intersecting particles is a central 
property of integral geometry |35| and follows from the 
IS0(3) invariance. 

Next, we transform {ei \ Cj*', 63*'') into the orthonormal 

coordinate system of the principal frame (vf' ,V2\n^^^), 
changing the notation for the normal direction 7? = 63 
to be consistent with Rosenfeld's and Wertheim's con- 
vention. The 3-dimensional cross product of the normal 
vectors 



ei = 7^1 = f2 A 7;3 



1 

S12 



77(2) X 77(1) 



(60) 



points now into the tangential direction of the intersec- 
tion. Combining the Euler form and the kinematic mea- 
sure, we obtain the intersection probability between two 
particles: 



1 1 



dDi A dD2 



riDixD2) 



rf I if (El n E2)d(Ei n S2)dL'2 (6i) 

™ Jv(D2) JSinS2 



47r 
1 



1 - C12 , ( 



471" Jv{D2) S12 



(2)1 

^13 + '^is ] 



dDo 



integrated over the intersection volume A = i^i n D2 and 
the kinematic measure with (j)i2 € T[D2). 

The transformation of the connection forms from the 
old reference system to the principal frame was done by 
Chern [36) . However, Wertheim's tensorial representa- 
tion [H] (see also [301 [Ml EOl [SS] ) has the advantage to 
be more closely related to Rosenfeld's definition of weight 
functions. In order to keep the discussion self-contained, 
we have included Wertheim's derivation in appendix [B| 
and present here only the result. 

Using the diagonal form of the Euclidean metric I and 
the curvature tensor K 



= vi ^ vi + V2 ^ V2 + n ^ n 

= Kl T^l (g) 7^1 + K2 V2®V2, 



(62) 



Rodrigues formula (33) yields the form: 

e\de-i = eiKei ds — ei — n <^ n) + A]ei ds (63) 
with the mean and tangential curvature 

K = ^(k1 + K2) , A = ^(ki - K2)(l'l ® i^l - i^2 (8) 1^2) ■ 

(64) 

With this change of notations and appendix [B| we fi- 
nally obtain Wertheim's representation of the kinematic 
measure 



1 



471" Jr(D2) JSiHEa 



ii:(Ei n E2)d(Ei n E2)dL'2 
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47r 



r(D2) J A=DinD2 



(I-n«(^n(2))(s(i)+;^(2)) 

(65) 



l + n(i)n(2) 



integrated over r^i G fi n D2 and r(L'2)- 

Now it is a simple task to expand the denominator in 
the geometric series 



1 



l + n(i)n(2) 



of tensor products and to rewrite the integral in the 
weight functions 



1 

47r 7r(_D2x...xZ5fc+i) 
xSinE2 



r(n2x...xDt+i) 

xStfc+i 



,(1),,(2) ,,(!),, (2) 



1), ,(2) V- 



(1) (2) 
-'AL+2'^(tL 



L=0 



+ (1 2) ^(3) ^ _ ^ ^{k+D ^3^^ ^ ^ ^ ^^^^^ 



(67) 



with the extended basis set of Rosenfeld's weight func- 
tions: 



ujal{D) 



1 
1 
1 



R{n)®^ 5{dD) 



A{nf^SidD) 



(68) 



47r 

uj,{D) = e(z?) . 
with the abbreviation: 

d{dD) = S{nr, dD) . 



(69) 



The normalization of the curvature dependent terms has 
been chosen to absorb the overall constant of Air. In the 
following we will see that these are all basis functions for 
3-dimensional, convex particles. 



4. The Three Particle Euler Form 

The third and last case is the Euler form for three in- 
tersecting particles. Its intersection Ei nE2 nSs consists 
of points, whose corresponding Euler class is a 0-form 
and independent of ojij . It therefore parallels the angu- 
lar dependent part of the Gauss-Bonnet equation ( 24 1 . 



As before (|53), the three normal vectors e^^\e^\e^^^ 



are converted into an orthonormal basis by the Gram- 
Schmidt method: 



Vl 
V3 



.(1) 
"3 

1 

M2 
VM3 



.(2) 



(er-(4^^«2).2-(4^^«r).i) 



(70) 



.(3) 



(3), 



and extended to the local frame 

Vl = Rij{li,l2,li)vj , e S0(3) , (71) 

interpolating between the three normal directions. Here, 
we can use the same argument that let to the simplifi- 
cation of ( 28 1 and replace the product of the connection 



forms by the volume of SO (3) in Euler angles: 

if(Ei n E2 n E3) = j sin (72)^71^72^73 • (72) 

However, 72 measures the angle between the vector and 
the a:2-axis and not the angle between the normal vectors. 
We therefore introduce a new coordinate system 



J2) 



' ^3 



J3) 



c(a3) 




(73) 



that is related to the Euler angles by 

c(72) = s{ai)c{a2)s{a^) +c(ai)c(a3) 
c(7i) = c{ai) , 0(73) = c(a3) . 



(74) 



The new representation of the Euler form ( 72 ) 

if(EinE2nE3) 

= j sin (ai) sin (02) sin {a^)daida2da^ 

= (1 - cos (012))(1 - cos (013))(1 - cos ((/123)) 

= (1 - n(i)n(2))(l _ r?(i)n(3))(l _ n^^)n(^)) 



(75) 



is a symmetric polynomial in the normal vectors. The 
remaining integration over the intersection space Ei H 
E2 n E3 reduces to a finite sum over its intersection points 



j ir(Ei n E2 n E3)d(Ei n E2 n Eg^ 

(1-Ci2)(l-Ci3)(l-C23) 



(76) 



St3 

X (5(n(iVA)(5(n(3)f^)5(^(3)f^) d^r^ 

1 1 - C12 1 - Ci3 1 - C23 

2 S12 Sl3 S23 ' 

pteSinS2nS3 23 
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where relation (B14) and the vector basis (73) for the 



normal directions has been used 



.S12S13S23 



(77) 



Furthermore, a factor 1/2 has been added to compensate 
for the double covering of the integration range, when 
instead of the Euler angles < 0, V' < Stt and < 9 < n 
the symmetric choice of the intersection angles 



< 012, 013, 023 < 27r . 



(78) 



is used. 

Next, we have to determine the kinematic measure 
with one of the three particles fixed in space. The deriva- 
tion parallels that of ( 59 ) and begins with the coordinate 
transformation of dD2 A dD^. Following the approach 
of [57], we rotate the locale frame of particle D3 by the 
matrix 



i?i(7i 




(79) 



in the 3 — > 1 direction and derive the new vielbein and 
connection forms for Di: 



,(3) 
^^13 

,(3) 
^^23 



,(1 



(1) 



C1WJ3 — S1UJI2 



-'23 



d-fi 



(80) 



0(3)^_,4i)+ei0W 



The same calculation has to be done for particle D2, 
where the matrix 



^23(72,73) 



generates a 2 — >■ 1 rotation 



C2 S2C3 S2S3' 
-S2 C2C3 C2S3 

-S3 C3 



(81) 



,(2) 



(1) 



23 = S2S30Jj_2 



(1) I (1) ^ 

5203^^3' + C2U}'23 ~ C2d73 



(2) (1) , (1) J 

^12 = C3^^12 + ^-i^lZ - "72 



(82) 



9(2) _ 
'2 



= -^2^1 



(1) 



C2C3(?W 



C2S3Bf 



The forms in the normal direction of D\ vanish by the 
constraint (32). We can therefore set the correspond- 
ing terms of ^g^'^^w^g^ and Wjs^ to zero and insert the 
transformed elements into dD2 AdD^. Performing an ad- 
ditional coordinate shift 72 -> 72 + 7r/2 and the change 
of basis ( 74 1 to transform from the Euler into the inter- 



section angles, we finally obtain the reduced kinematic 
measure 



dD2 A dDs 

= (si2Sl3S23)^C?Sidi;2rfS3d0i2d0i3d023 



(83) 



with the kinematic measure of the surface dE defined in 



(40). 



Collecting terms, the Euler form ( 76 ) intersecting with 



k — 2 further particles is determined by: 
1 

in 



r{D2X...xDk + i) 



Ki^i n E2 n S3)d^2 . . . dDk+i 



{pt} 



1 

8^ /r(D2X...x_Dfc + i) 

xStfc+1 



X (1 - n(i)n(2)) (1 _ ri(i)n(3)) (1 _ n^^)n^^^) 



(84) 



r(2) 



(3), 



and can be rewritten in the basis of the weight functions, 
defined in ( 68 ) , after expanding the product of ( 84 ) : 



1 

An 



r{D2X...xDk + i) 



Ki^i n S2 n S3)rf^2 . . . dDk+i 



{pt} 



\, /I), /2) (3) 
^ ^aO ^aO ^crO 

k + l) I 



1 

87r jT{D2X...xD, 
xStfc+i 

(1) (2) (3) (2) (1) (3) (3) (1) (2) 
, ,(1), ,(2), ,(3) (2) (1) (3) (3) (1) (2) 



,,(1),,(2), (3) 
^a2 ^cr2 ^a2 



J^K..J^+^U\Ad^2-.-d^k+l 



(85) 



As required, the result is invariant under cyclic permu- 
tations of the indices (1, 2, 3). 



For the first two integrals (67 47) it was possible to 



scale the pre-factor to one by a suitable definition of the 
weight functions. The same is not possible for (85), as 



it depends only on the previously defined weights. The 
three particle integral has therefore an overall pre-factor 
of l/87r. 



The three intersection probabilities (67 47 85) are 



complicated polynomials in the weight functions. How- 
ever, here we have shown, by explicit calculation, that 
these three cases are all we have to consider under the 
given restrictions on the manifolds. The five different 
types of weight functions ( 68 ) are complete in this sense 



and provides the basis for higher loop orders. The group- 
ing of the weight functions into five classes can be stated 
more formally by their scaling dimension under the co- 
ordinate transformation r — >■ Xr. 

Let us summarize the results of this section: 

Theorem II. 1 The Euler form lo^ of the kinematic mea- 
sure of a stack Stk of 3-dimensional, convex Riemannian 
manifolds decomposes into a symmetric sum of weight 
functions 

uj^i^i n . . . n Sfe) = CAi...Ak ujAii^i) ■ . . WAfc(Sfe) 

io^iJ:in...nDk-inDk) (86) 
= uj^{Y.in...r\Dk-i)LjyiDk) 
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where an implicit summation over the multi-index A € 
{x, V, kL, AL} for L = 0, 1, 2, . . . is understood. The nu- 
merical values of the coefficients CaiA2... follow from { 61 



They depend on the dimension of the embedding 
space and the particle but are otherwise independent of 
the manifold's geometry. 



The weight functions ( 68 ) provide a complete basis set, 
in which the intersection integrals can be expanded. They 
are unique with respect to the Euler form. Their scal- 
ing dimensions group the weight functions into four sub- 
classes: 



[u]} = 3 , [ul^] ^ [c^Xl] - 2 , [ul^] = 1 , K] = 



(87) 



III. 



RESUMMATION AND THE ROSENFELD 
FUNCTIONAL 



A. The Functional of Rosenfeld and Tarazona 

1. Rosenf eld's Three Postulates 

The local decomposition of the kinematic formula for 
one, two, and three particle intersections clarifies the 
mathematical aspects of Rosenfeld's approach. However, 
it remains to combine the resulting weight functions into 
the free-energy functional. A first naive attempt of in- 
serting the reduced virial integrals into the corresponding 
expansion of the chemical potential 



n=l 

PJ"^ J n{r)dp{r}d^r 

fails. The reason lies in the decoupling of the particle 
density p from its geometric properties uja that allows to 
add a particle by the integration of ( 88 1 without adding 
the particle's volume To find a corresponding gener- 
alization, let us reconsider Rosenfeld's derivation of the 
fimctional [T^ (see also [3]). 

The infinite number of weight functions (|68|) reduces 



to a finite subset for spheres, whose principal curvatures 
Ki = Ki causes the A-dependent terms to vanish. The 
second virial integral of a mixture of hard spheres with 
M components reduces therefore to a finite sum of only 
six weight functions. 

Ai,A2 



[i ^ i) 



(89) 



tensor product is a short form of the convolute integral 



^M^^A,^ I i^M{^A-ri)uJ^AA^A~rj)d^rA (90) 
JDinDj 

depending on the particle positions fi , fj in the embed- 
ding space M.^ and the intersection point rA ^ Di O Dj. 
From the decoupling of the integral measure ( 89 1 into 



single particle contributions follows the splitting of the 
entire second virial integral, weighted by the 1-particle 
densities Pi{ri): 



/3i(A,i^,) 



J2 C^^^^ / p(*)p(j)K, ®Orf7.rf7jrf'''A 
C^'""' I nAArA)nA2irA)d\A 



written in the weight densities: 



M . 

nA{rA) = ^ 



Pi{fi)u\{rA ~ n) d-ji 



(91) 



(92) 



As has been discussed |II A[ the pairing of one weight 
function with the 1-particle density is a consequence of 
the single intersection domain of the second virial clus- 
ter. However, it is natural to generalize this construc- 
tion further to particles with k intersection centers. The 
corresponding integral then combines k weight functions 
with the 1-particle density: 



"■Ai,...,Afc(r4i,..-,r4j 



M . k 

y2 Ptin)Y\ ^A^ i^A,. ~ ^0 
7^1 J no.) f:=i 



(93) 



generalizing the 2-point densities of the exact third virial 
integral Q. Such "fc-point densities" are the central ob- 
jects in analyzing higher loop diagrams. With increasing 
loop order increases also the order of the /c-point densi- 
ties. This can be seen by assuming that all g loops begin 
and end at the same particle. The loop diagrams then de- 
couple into sets of fc-point densities for 2 < fc < 2g. The 
only diagrams that contain 1-point densities are therefore 
the intersection stacks of 17 = 0, as has been explained 

From the observation that the leading contribution of 
the free-energy factorizes into products of weight densi- 
ties, Rosenfeld postulates three assumptions about the 
structure of the functional: Firstly, the free-energy is an 
analytic function in the weight densities, i.e. it allows a 
polynomial expansion in ua 



where i, j = 1, . . . , M runs over all types of spheres. The 



/?J-(M)= / <^^^{[nA])d\ 



(94) 
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Of course, we have seen in section II A that this assump- functional (97) by this term alone is in excellent agree- 



tion is not true in general. However, the functional form 
of can be further restricted by observing that the in- 
tegral (94) has to be invariant under coordinate scaling. 



The second assumption is therefore that the free-energy 
functional is a homogenous polynomial under the trans- 
formation r — > A^^r with the scaling dimension 



(95) 



of the free-energy. The possible combinations of weight 
functions are therefore constrained by their scaling di- 



mensions ( 87 1 with the exception of the scale indepen- 
dent LUy-. 



+ fi{nv)nl„ + f5iny)n^on^ina 



(96) 



With the third postulate, Rosenfeld further assumes that 
the functional is a solution of the scaled particle differ- 
ential equation [3l [11] . In this way it is possible to deter- 
mine the dependence of the unknown functions /i , . . . , /s 
on the scale-invariant weight density Uy . The free-energy 
functional is then known up to the integration constants 
of the solutions of the differential equation. For /i, /2, /s, 
they can be read off from the second virial contribution; 
but the constants for /4 and have to be determined 
by comparison with analytical results obtained by alter- 
native methods. The functional has thus the preliminary 
form IMl: 



*preiim(K]) = In (1 - riy) + 



1 -Uy 



+ 



1 r4o - 3ncronain„i 



247r 



(l-n,)2 



(97) 



Later on, it has been shown that this functional leads 
to an unphysical singularity, when the positions of the 
spheres were constrained to lower dimensions [121 HH]- 
The source for the occurring divergence is the third term 
in the functional. This led Rosenfeld and Tarazona to 
look for alternative third order polynomials compen- 
sating the singularity. Several suggestions were made 
[TOl fT71 and compared to simulations. The most 
promising modification today is Tarazona's [17^ replace- 
ment: 



1 

16^ 
1 

16^ 



(98) 



;(l-Ci2)(l-Ci3)(l-C23)-A/3 



with Afs from (51 ). Comparing this semi-heuristic result 
to equation (75), identifies the first term as the three- 



particle intersection probability of the stack. In [121 HI] it 
has been shown that the corresponding correction of the 



ment with simulation data of the bulk-fluid free-energy of 
hard spheres. The fluid phase is therefore well described 
by the intersection probability of stacks. However, it has 
been shown in [16, that the Lindcmann ratio for the fcc- 
lattice is underestimated by this functional. This is cor- 
rected by the second part of ( 98 1 , improving the equation 
of state for the solid region 17J. In the next section we 
will argue that this term is part of the 1-loop correction 
of the third virial diagram. 

The final form of the Rosenfeld functional for hard 
spheres pT] is obtained by replacing the third term of 
( 97 ) by Tarazona's expression ( 98 ) : 



*R (Kl) = In (1 - Tly) + 



"•ctO"-ct1«c71 - n<jirtCT2"-CTl + ncr2«<T2"o-2 ^ »^(T0"cr2«<T2 



il-nyf 



(99) 



This result provides one of the currently best approx- 
imations of the fiuid phase structure of hard spheres, 
only surpassed by the White Bear version [211 113] . How- 
ever, this improvement has been obtained by adjusting 
the functional to simulation data, whereas the correction 
( [98| is geometrically motivated. Apart from the Afs-term 
in (|98[), we have already derived all of its contributions 
and pre-factors from the 0-loop order. 



2. Replacing the Scaled Particle Differential Equation 

The chemical potential enters the fundamental mea- 
sure theory via the scaled particle differential equation. 
Its origin is a semi-heuristic relation between the chemi- 
cal potential and the pressure /z™ — pvi in the low den- 
sity limit that becomes exact at diverging particle volume 

— >■ oo. This limit allows to relate the chemical poten- 
tial of the free-energy T to the pressure representation 
of the grand potential —pV — VL ^ T — piST /5pi. Intro- 
ducing the functional derivative: 



Spijn) 



Sij S{fi - r-j) 



(100) 



which selects the weight function when applied to a 
weight density 



Spj{rj) 



riA 



{ta) = j ^^^a{'^a~ ri)Sij S{ri - fj)tf 



uJ^Ai^A-rj) , 



(101) 



the chemical potential nf^ of the free-energy functional 
has the form: 



/3Mr(^r;) 



E 



^Piiri) ^ driA Sp,{fi 
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dUy 



9$ 



= {-^W + T.''A^+n^)u;Ur,r^) , (102) 



The two derivatives with respect to pi and ua = Pi^\ 
are of course not independent from each other and do 
not commute ^ p.yi. It is therefore important not to 
interchange the order in the integration 



(107) 



J^"""^ = / \iiySny = / {fi^y 6pi)Sny 



assuming that all contributions of lu^ vanish in the 
Vi ^ oo limit except for w,^ . From this follows the scaled 
particle differential equation: 



R 



dUy 



E 

A 



HA 



duA 



(103) 



The arguments leading to this result are by no means 
trivial: The scaled particle limit allows the identification 
of the particle volume Vi as the embedding volume V, 
resulting in the unpaired index v in the last two lines of 
(102 1. Another striking feature is the dependence of the 



chemical potential on the two different coordinate sys- 
tems of the particles fi e Di and those of the intersec- 
tion region r € Stk ■ This indicates a further difficulty in 
identifying the chemical potential as an external potential 
coupled to the particle density. To obtain a symmetric 
formulation in the densities pi and ha, let us define the 
chemical potential for the particle volume n„: 



^yir) :=/3 



S:F-^{[nA]) 
5ny{r) 



(104) 



In principle it is possible to define an infinite set of 
chemical potentials for the weight functions uj\. How- 
ever, "^y is the only physically relevant one. This can 
be realized in two different ways: Firstly, 5ny is again 
scale invariant, which follows from [p,] — —[(Pr] — i and 
[uy] = 0. '^y has therefore the same scale dependence as 
the free-energy. This complies with the interpretation as 
the energy change by inserting a particle into the system 
and the observation that ujy is the only scale invariant 
weight function. Secondly, it follows from (45) that the 
intersection probability of a stack St^ of order fc > 3 will 
only change by a factor w*, when an additional particle 
is inserted. This corresponds to a formal integration over 
LOy coupled to the particle density pi. 

The functional derivative (104) can be inverted by in- 
tegration 



(105) 



^j7cx ^ / ^^^^Jf)Sny{r) := / ^yduycfr 



= j $'=='(r)dV 



and relates the chemical potential to Rosenfeld's free- 
energy density. It also allows a natural interpretation of 
'^y as the integral of the functional derivative 



S T^^Hua]) 
5pi{n) Sny{r) 



(106) 



Now, piy has the right structure for generalizing the 
virial expansion ( 88 ) to the weight function depending 
terms /3„(lli^)p". Furthermore, it is extensible to ar- 
bitrary loop orders. Inserting the expansion (88) into 



( 107 1 with subsequent integration over pi gives a general 



relation between the virial expansion and the free-energy 
density ^10^■. 



oo 

$-(M,f) = cn.+^— / 
fc=i •' 



p'^+^Pkduy . (108) 



The integration constant c is itself a functional of the 
remaining weight densities ua for A ^ v to be determined 
by comparing to the low-density limit. However, 
the scaling dimension restricts the possible dependence 
to c cx n^, with a universal constant to be determined in 
the next section. 



Equation (108) generalizes the virial expansion (88) of 



the free-energy to the functional form depending on the 
weight densities. It is an exact relation and independent 
of the semi-heuristic scaled particle theory. Once the 
virial coefficients are known, we can derive the functional 
by a simple integration over ?i„ for any loop order. 



B. The 0-Loop Order of the Free- Energy 
Functional 



With the derivation of the intersection probability of 
particle stacks ( 86 1 and the virial expansion of the free- 
energy in terms of the weight densities (108), we can 



finally put the pieces together and prove our hypothesis 
(12) that Rosenfeld's functional <f>^ corresponds to the 
leading order $o of the loop expansion ( 14 1 . This is done 



in two steps: deriving the virial integrals for any diagram 
of zero order, and then adding them up into a generating 
function. 

In section |II A| we have seen that a Mayer cluster of 
loop order g decomposes into a series of topological dia- 
grams 



(109) 



n=0 



of which the leading order corresponds to the intersec- 
tion probability of a stack St^+i. Following the discus- 
sion from section [II C[ the corresponding cluster integral 



1 cr 



/l,2 • ■ • fk,k+l dji . . . djk+1 



r(_Di 



16 



kl Jr{D2X...xDk + i) 



(110) 



is identical to the averaged Euler form, integrated over 
the kinematic measure of fc + 1 particles. Here we have 
used that the symmetry coefficient is cr = 1 and that the 
volume factor V cancels after integrating over the coordi- 
nates of the center of gravity. In principle it is possible to 
extend the integral to mixtures of particles by including 
an additional index. However, this is not necessary, as 
the final result will depend on the weight densities ( 92 1 , 



which automatically include the right combinatorial fac- 
tors. We can therefore restrict the discussion to a single 
class of particles without loss of generality. 

The boundary of a stack of identical, 3-dimensional 



particles has been derived in ( 44 ) and reduces to the sum 
of three contributions. The branching rules of (861 can 



then be used to algebraically split the Euler form of ( 110 ) 
into the volume dependent weight functions 



+ fc(fc + i)wx(sni;)w^~i 

+ k{k + i)(fc - 1) uj^{^ n s n s;)a;^'^^ 



(111) 



and further into the decoupled product of weight densi- 
ties: 



UJ^CE n S n S) = CAiA2A3i^Ai^A2^A3 

a;x(S n S) = CaiA2 i^Aii^Aa 



(112) 



where an implicit sum over the paired indices is under- 
stood. We also introduced the trivial constant = 1 
to keep the notation symmetrical. In anticipation of the 
following derivation of the Rosenfeld functional (99), it 



is useful to separate the dependence on the highest and 
lowest weight functions , ujy from the Euler form and 
to introduce the index notation 



A = (x, V, a) = (x, w, kL, AL) 
deduced from Theorem III. II 



(113) 



integral for a stack 



Inserting (fill and (112) into (110) yields the virial 



= + 1) / \Cxv ^x^v +kCa^a2^ai^a2^'; 



fe-1 



fc+1 



+ k{k - 1) c„ 



4=2 



(114) 



of fc -I- 1 indistinguishable particles. The virial coefficient 
is a homogeneous polynomial of order fc -t- 1 in the weight 
functions and combines with the particle density p^^^ to 
a polynomial of weight densities. Inserted into (108), we 



obtain the result: 



<{[nA\) = c + Y.^/'~'Pl 



k=l 
1 



, 1^ \ — 1 1 4- ^ ^ai ^a2 



(115) 



2 a 



(1 - n„) 



3 ■ 



The integration constant c can now be uniquely deter- 
mined by comparing it to the ideal gas limit, where the n„ 
dependence has to vanish. Inserting the value c = C^v^x 
and integrating over the n„ density gives the final excess 
free-energy functional of the 0-loop order: 



^"^{[nA]) 



-C- 



,y In (1 — Uy) -\- — 



(116) 



^Ql 

■ (l-n„)2 



Comparing this result to the Rosenfeld functional (97 1, 



we have finally proved our hypothesis (12). 

This result also allows a formal extension to 



D- 



dimensional particles embedded into the odd dimensional 
M.^ . Because the Mayer expansion is independent of the 
dimension of the physical system, nothing will change by 
this generalization. Extending the boundary stack (44) 



to D dimensions and the corresponding splitting of the 
Euler form ( 86 ) results in a free-energy functional 



D 

E 



...ctk ^^ai 



gk 



(117) 



(n.u) = (1 - n„) In (1 - n„) + n„ 



that can conveniently be written by the generating func- 
tional 4>. The same observation has been made before in 
|20) . where 4>{ny) has been derived in the freezing limit, 
when the particles are located in caverns. Here, we can 
see that the generating functional carries the volume de- 
pendent parts of the boundary of the universal stack as 
defined in ([T]) . The Rosenfeld functional has now the sim- 
ple interpretation as the intersection probability of USt. 

Thus we have shown that the 0-loop order of the virial 
expansion leads to the Rosenfeld functional. However, it 
only reproduces the first term of Tarazona's correction 
([98]). Therefore, one might guess that the Ms-dependent 
part belongs to the 1-loop correction of the third virial 
order (|4| as will be investigated in a subsequent article. 



IV. DISCUSSION AND CONCLUSION 

In this article it has been shown that the Euler form 
K(dStk) determines the intersection probability of a par- 
ticle stack of order fc and that its generating function 
reproduces Rosenfeld's functional. These results explain 
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and generalize Rosenfeld's previously unproven observa- 
tion [lTl[39l|40j that the second virial integrand is related 
to the Gauss-Bonnet equation. For two intersecting con- 
vex particles the results of Werthcim [IT and Hansen- 
Goos and Mecke [301 HB] are confirmed by cxplicitely de- 
riving the Euler form from first principles. However, go- 
ing beyond the second virial, we further derived the pre- 
viously unknown Euler forms for fc > 3 and their splitting 
into weight functions. 

Motivated by the success of Rosenfeld's functional for 
the liquid region, we made the Euler form the foundation 
of the fundamental measure theory and its extension be- 
yond the currently known functional. It has been shown 
that the Mayer clusters of hard particles split into inter- 
section diagrams that can be classified by their number 
of loops and intersection points, where the latter corre- 
sponds to a particle stack. The leading contribution, the 
0-loop order, is then the only part of the free-energy that 
can be represented by a functional with only one inter- 
section point. 

From this follows that the fundamental measure theory 
allows the systematic derivation of the free-energy func- 
tional for each loop order; a result that is in fundamental 
contrast to DFT in quantum mechanics, where the devel- 
opment of a functional is only restricted by the existence 
theorem of Hohenberg and Kohn [SH]. This property of 
hard particle physics is probably a consequence of the 
invariance of the Euler form under geometric deforma- 
tions. As long as the homotopy type and therefore the 
topology does not change, we obtain the same functional 
form. And even if we include complex geometries like tori 
or hollow spheres, the additional terms still derive from 
an Euler form. The only constraints we have to consider 
are of physical nature and are related to concave geome- 
tries. 

The infinite number L — 0,1,... of tensorial weight 
functions provide a practical problem in the calculation 
of higher loop orders. Since we cannot derive an infinite 
set of integrals, it is necessary to stop at a certain or- 
der. A first hint gives Wertheim's calculation of the third 
virial integral for prolate and oblate spheroids [131 HI] . 
He shows that the aspect ratio A < 10 differs from the 
simulated result by less than 3%, when the L < 2 terms 
are included. This indicates that the expansion of the 
denominator 1 + n^^-'n^^^ is fast converging for most of 
the physically interesting cases. 

Also of importance is the influence of the number of 
loops and intersection points. As explained in section 
|II A[ each intersection point of a diagram is dressed by 
the universal stack, as shown in FIG. [2] whose free- 
energy contribution is already known from the 0-loop 
order. Consequently, each intersection carries a factor 
of (1 — and (1 — n^)~^. From this follows that the 

divergence of the resummed third virial integral *I'(i,3) 
of FIG. [5] is at least of order (1 — n„)~'^. The influence 
of diagrams decreases therefore significantly with their 
number of intersection points. We therefore expect no 
new physical effects by including higher intersection or- 



ders. This is consistent with our hypothesis that only 
higher loop orders correspond to long range effects be- 
tween particles, as indicated by the generating function 
of all 1-loop diagrams. 

Another aspect worth considering is the dimensional 
influence of the particles and their embedding space. If 
the codimension is larger than 1, the particles do not 
necessarily intersect, while approaching each other. The 
mathematical formulation is then more complicated and 
requires the introduction of equivariant differential forms 
|54| : in the physical literature this is known from BRST 
quantization [50] . We have also seen that the Euler form 
vanishes for odd dimensions and gets replaced by higher 
order invariant forms. This is a consequence of the Bott 
periodicity [BT and offers a direct link between the math- 
ematical and physical properties. It is even possible that 
this relation can be further extended to a more detailed 
understanding of the relation between topology, geom- 
etry and the physical phase structure of particles. For 
example, one might ask, if the geometry of a particle and 
its mixtures can be tested by their phase diagrams? 

An important step in this direction is the numerical 
calculation of weight functions and the minimization of 
the grand potential functional For the 3-dimensional 
particles it is possible to reduce the problem to a tri- 
angulation of the surface and to replace the connection 
form by a sum over the outward angles, analogously to 
the derivation of the Gauss-Bonnet equation. The result- 
ing polyhedrons are then placed into a Voronoi diagram, 
whose boundaries are varied until the minimum of the 
free-energy has been obtained. This approach would al- 
low the analysis of even more complicated particle dis- 
tributions than the isotropic or periodic structures inves- 
tigated so far. In addition, it would also allow a better 
understanding of the origin of phase transitions. For in- 
stance, the particles in the nematic and smectic phase 
are parallel oriented, minimizing the 0-loop contribution 
of the free-energy by setting one or more of the inter- 
section angles to zero. However, understanding such ef- 
fects requires the derivation of higher loop orders and will 
therefore be postponed to the next article. 
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Appendix A 

It is enlightening to compare the local formulation of 
Chern to the approach of Minkowski ^SSI , which was 
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the basis for the calculation of Isihara and Kihara [51 [5] ■ 
We will therefore give a short summary of their deriva- 
tion that led to the first general equation of the second 
virial coefficient of convex particles. Let pi G Di be the 
coordinate vector of the two convex particles i = 1,2. 
The excluded volume under translation and rotation of 
the particles is then calculated by first deriving the dif- 
ferential volume element dVi2 of the shifted coordinates 
followed by the rotational averaging. We first obtain 

1 



dVi2 = j^d^{pi+P2) 



[d^pi + Mp2 A d^pi + (1 O 2)) 



(Al) 



= dVi + -d[p2,dpi,dpi] + (1 o 2) 
= dVi + H2dSi + (1 o 2) 



with an implicit integration in the second part and the 
support function H = pe^. The orientation has been cho- 
sen such that the normal surface vector of particle D2 
at contact is —63. This allows to simplify the determi- 
nant, indicated by the square brackets, via the relation 
[p2, dpi, dpi] = [p, 6'"eQ,^'^e^] = 9i f\ 92{pez) as shown 
in |46]. The rotational averaging over the coset space 
SO(3)/SO(2) reduces again to the multiplication by the 
connection form uj^ A LO2 — kiK20i A 02 — KdS: 



J {dVi2),ot = J K2dS2 A dVi 



J H2K2dS2 A dSi + (1 ^ 2) 



(A2) 



The product between the support function and the Gauss 
curvature can further be simplified by the substitution 

m 



= y d[p, 63, ^63] = J [dp, 63, ^63] - [p, des, des 
= 2 J [HK - M)dS . 



(A3) 



Inserting into equation ( A2 ) , finally gives the result of Isi- 



hara and Kihara as a special case of Minkowski's formula 

ESI 



1 

47r 



{dVi2),,, = X2V1 + ^7*251 + (1 O 2) . (A4) 



This result can also be obtained in a coordinate-free rep- 
resentation by the Lie-transport exp {Cx2)dVi — dVi2 of 
the volume form and Stokes formula 



Cx,n2 = / d{ix,n2) = / ^Xl^^2 . (A5) 

D JD JdD 



Appendix B 

In the following, we will give a short account of how 
to transform the two particle Euler form ( 56 1 to the co- 



ordinate dependent representation ( 65 1 of Wertheim, as 
used in ^41j . 



The Euclidean metric ( 16 ) in the orthonormal principal 



frame {iyi,i'2,n) is the diagonal tensor 

Vij = ei ® e.j = lij 

= {vi®vi + V2®V2+rt® n)i 



(Bl) 



of ( 62 1 . The related connection tensor ( 62 ) then follows 
from the exterior derivative of the normal vector 63 = n: 



de^ — idSaea = Ka9a ® = l^a&a ' 
— {niVl ® Vl + K2t^2 ® 1^2)^ ds 

= K ids 



I Cadp 



(B2) 



using Rodrigues formula ( |33| , the representation of the 
vielbein 6a — e-adp, and by observing that the tangential 
vector at each point p G Si H S2 lies in the direction of 
t ^ n'^) X ■nP'\ The derivative dp — tds therefore is the 
differential line element ds pointing into the direction of 
t. 

In order to separate the normal vectors from the prin- 
cipal frame, Wertheim rewrites the connection form |41j : 



(kii7i (g) i7i + K2i'2 ® 1^2) 



(I 



n ^ n — V2 ^ V2) 



-K2[l — nii$ n 



- 1^1 vi) 



(B3) 



= -(ki + K2)(II -n®n) 

+ ^ (k1 - K2){vi ®Vl-V2® V2) 



with the mean and tangential curvatures defined in ( 64 1 
The connection then yields the form 



'^is = eides = tKtds 



(B4) 



of (56 1. In a second step, the normal vector n*^^^ is sepa- 



rated from the curvature depending parts of particle 1: 

X K(i) (nW xn(2)) 



= -n*^^ X n''^^ {kiVi ^ vi + K2i'2 ® V2)n^^^ x r?^^^ 



x vi ® vi x n 



(1) 



= n^^^ {ki1^2 (X" 1^2 + H2'^i (gi i'i)n^'^^ 
= kJ^) n(2) , 



(B5) 



using the orthonormal relation vi x 1/2 — n and introduc- 
ing the adjoint connection tensor: 



= K (I ■ 



n ^ n) 



(B6) 
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Inserting these results into ( 56 1 

1 - Ci2 (1) _ 1 - Ci2 



13 



Sl2 
1 - Ci2 





Sl2 


1 


- Ci2 




s2 
■'12 




1 


1 


+ Cl2 




1 


1 


+ Cl2 




1 


1 


+ Cl2 



Sl2 
Sl2 



^(1)- 



Sl2 



Si2 



Sl2 



r(2)r^(i) 



7(2) 



ds 

Sl2 



;'^^'Hl-c?2)-n(2)A(i)ri(2); 



Sl2 



^(2)^(1)^(2) 



1 + n(l)r?(2) 



ds 



\nW X n(2)| 



and using the integral representation by (5-functions 



1^131, ,(1) 



Sl2 



'13 



DinD2 



n(2)A(i)r?(2) 



(B7) 



l + n(i)n(2) 

this reproduces the first part of Wertheim's equation 



(65 1. The second part follows accordingly by replacing 



the particle indices 1 o 2. 

The integral representation used in (B7| extends the 
integration along the line element ds to the entire embed- 
ding space. This and similar relations are readily derived 
from the linear coordinate transformation 



1] = np , C = mp , ^ = eix + e2y + e^z 



(B8) 



at the point p — (x,y,z) and its corresponding Jacobi 
determinant: 



drj A d^ A d^ — \ det (n, rn,e} \ dx A dy A dz 
— \n X m\ d^p . 



Applied for the integral of an arbitrary test function F 
and two (5-functions 



F{pj 5{np) S{mp) d^p 



d-q dC d£, 
In X ml 



(BIO) 



d^ 



it reduces to the line integral along ^, as used in equation 
(iBTl). 



With one (5-function included, the corresponding trans- 
formation 



T] = np, ( = = eix + e2y + e^z 
and e A e = e yields the result: 



(Bll) 



F{p) 5{np) d^p 



drj dC d£, 

\en\ 

F{C,0 dS 



(B12) 



with det (n, e , e) = en and the differential surface ele- 
ment dS„ in the outward pointing n direction. 

Analogously, the integral of three (5-functions reduces 
to a sum of intersection points {pt} in the variables 



r] — rip , C — mp , ^ — Ip 



(B13) 



solving the algebraic equation ?/ = C = C = 

/ Fip) 5{np) 5{mp) 5{lp) d^p = J2 ,J^^^lr, (^^4) 
J ^p^j|(nxm)Z| 

as appears in the equation of the intersection probability 



(B9) of three particles (761 
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